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Abstract 

We study the effect of electron-electron (e-e) interactions on compressibility, 
capacitance and inverse compressibility of electrons in a quantum dot or a 
small metallic grain. The calculation is performed in the random-phase ap- 
proximation. As expected, the ensemble-averaged compressibility and capac- 
itance decreases as a function of the interaction strength, while the mean in- 
verse compressibility increases. Fluctuations of the compressibility are found 
to be strongly suppressed by the e-e interactions. Fluctuations of the capaci- 
tance and inverse compressibility also turn out to be much smaller than their 
averages. The analytical calculations are compared with the results of a nu- 
merical calculation for the inverse compressibility of a disordered tight-binding 
model. Excellent agreement for weak interaction values is found. Implications 

for the interpretation of current experimental data are discussed. 
PACS numbers: 71.55.Jv,73.20.Dx,71.27.+a 

Typeset using REVTgX 



I. INTRODUCTION 



The interplay between disorder and interactions is one of the most exciting topics in meso- 
scopic physics. Both factors play an important role in determining the change of the number 
of electrons in a mesoscopic system as the chemical potential of the system is changed. The 
investigation of this problem is especially timely since ground-state properties of disordered 
interacting systems have recently become experimentally accessible. For example, by mea- 
suring the spacings between the gate voltages for which a quantum dot connected to external 
leads exhibits conductance peaks, one can infer the inverse compressibility of the electron 
gas in the dot Devices of this type have recently been used to measure the many- 

particle ground-state energy as a function of the number of electrons in a dot and large 
fluctuations (compared to the single electron level spacing) in this property were observed 
|3|-[5|]. Large fluctuations have also been observed in earlier experiments in which the inverse 
compressibility of an insulating indium-oxide wire was measured || . 

In the absence of interactions the derivative of the chemical potential fi with respect 
to the number of electrons iV (i.e., the inverse compressibility, equal also to the second 
derivative with respect to iV of the ground state energy) of a system equals to the single- 
electron level spacing. It is well known that the single-electron level statistics in disordered 
systems are connected to the statistical properties of random matrices [0-0- The level 
spacings exhibit fluctuations and therefore one expects that the inverse compressibility of 
different samples, as well as the inverse compressibility for the same sample at different 
chemical potentials, will also exhibit fluctuations of order of the single-electron mean level 
spacing A. In the non-interacting regime there are many theoretical approaches which can 
be used to calculate these fluctuations, among them are the random matrix theory (RMT) 
@, perturbative diagrammatic calculations || and the supersymmetry method fliO| , [TT]| . 

On the other hand, once electron-electron (e-e) interactions are included the situation 
becomes more complicated. In contrast to the high energy excitations of an interacting 
system for which a RMT description works very well (for example the high excitations of 
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a nuclei || and of disordered interacting systems [12j ), The ground-state energy of an 
interacting system as a function of the number of particles can not be understood by means 
of RMT. 

The average compressibility (dN/dfi), as well as the average capacitance (C), are ex- 
pected to decrease as result of interactions since it costs more energy to add particles into the 
system. This is somewhat similar to the situation for the average polarizability suppression 



by the e-e interactions flqjli 



It is not a priori clear how the interactions will influence the fluctuations. Previously 
we have calculated the influence of e-e interactions on the sample-to-sample fluctuations 
of the polarizability ||15|| . By analogy one might expect that the relative magnitude of the 
compressibility and capacitance fluctuations should be suppressed. 

The average inverse compressibility (dfi/dN) is expected to increase as a function of 
interaction strength for repulsive interactions. From the experimental evidence it seems 
that the fluctuations in the inverse compressibility is proportional to (dfi/dN) 0-§J. From 
numerical calculations one learns that for strong interactions the mean root square of the 
inverse compressibility are proportional to its average ||, with a proportionality constant of 
about 0.15 for a very wide range of interactions and disorder strengths. However, it is evident 
that for weak interactions this ratio can not be constant. Our goal is to study the behavior 
of the inverse compressibility for weak interactions, and to understand the transition to the 
strong interaction regime. 

In this paper we extend the perturbative diagrammatic calculations to include e-e inter- 
actions. The average compressibility, capacitance and inverse compressibility as a function 
of the interaction strength are calculated. The fluctuations in the compressibility and ca- 
pacitance are found to be strongly suppressed by the e-e interactions. The fluctuations in 
the the inverse compressibility of the sample can also be deduced from the fluctuations in 
the compressibility. It turns out that the fluctuations in the inverse compressibility, while 
not suppressed are relatively small. 

This is checked by an exact diagonalization numerical study of an interacting, spinless, 



disordered tight-binding model for nearest-neighbor interactions as well as for long-range 
Coulomb interactions. For weak interactions the numerical results follow the perturbative 
calculation. As may be expected from earlier studies, || deviations appear as the strength of 
interaction is increased. Thus, a transition from weak-interaction behavior of the fluctuations 
to strong- interaction behavior is evident. In the end we discuss the significance of these 
results for the interpretation of the experiments. 



II. AVERAGED COMPRESSIBILITY 

The compressibility may be written in terms of the exact Green function of the system 

as: 

f. = -I » to TV G*( £ ) (1) 

where /i is the Fermi energy of the system, N is the number of particles. Given the Hamil- 
tonian H, the exact Green function can be written as 



For the non-interacting case Eqs. (|l|J^) take the form 



dN 1 d T . / 

— — = — — Im / de dr& (r,r,e), (3) 



<9/i 7T dfl 

where Q is the volume and 

G R (f,f',s) = (AG\s)\r>) = ^: { ^2 ' ( 4 ) 

Here ip n {f) (e n ) is the n-th eigenvector (eigenvalue) of the system. This defines the density 
of states at the Fermi energy v 

According to Eq. ([Sp in order to add one electron to the system the Fermi energy should 
change by A. This can be also represented in a diagrammatic perturbation theory by the 
diagram shown in Fig. |l](a) corresponding to 



/dN\ = %_ r» r d ^ G R,^ F) £ )G R (f, f, e)) - (G A (f, f, e)G A (f, f, e)) (6) 
\ a/i / 7r j-oo in 

where (...) denotes an average over different realizations of disorder. 

Now let us discuss the influence of e-e interactions on the averaged compressibility. One 
would expect that repulsive e-e interactions will reduce the compressibility since it costs 
additional interaction energy to insert an electron into the system. We begin by investigating 
the role of short range interactions represented by 

U(f, r') = Xu-^if- r') = a d U5(r- f ') (7) 

and 

U{q) = \u- 1 = XQA = a d U, (8) 

where A = e 2 z/a d_1 is a dimensionless coupling constant and U = e 2 /a; here a is the range of 
the interaction and d is the systems dimensionality. The compressibility in the interacting 
case may be represented by the diagrams shown in Fig. |l|(b), corresponding to 

where xio) ^ s the result of summing the diagrammatic series shown in Fig. |l](c). For a short 
range interaction 

= TTT (10) 

Inserting the result for x(o) m t° Eq. ([5]) one obtains 

(ijl) = (i + a) (a) = (i + a d u/nA) (a) • 

Thus, the additional energy needed to add an electron due to the e-e interactions (the 
charging energy) is a d U/Q. This is exactly the result expected if one assumes a constant 
density of electrons in the system. 

The calculation of the compressibility for the long-ranged Coulomb interaction can be 
performed along similar lines. For the Coulomb interaction 



and for an infinite system 

U{q) = S d e 2 /q d -\ (13) 

(where §2 = 2tt and S3 = 4ir) which for future convenience in the numerical simulations 
may also be written as U(r,r ') = Ua/\f — f '\. As can be seen in Eq. (||) the zero mode 
(q = 0) determines the behavior of the average compressibility. Therefore, one must treat 
the zero mode carefully, since simply inserting the infinite system value of the interaction 
for q = will lead to no compressibility, i.e., infinite capacity. The zero component of the 
Fourier transform for a finite system is equal to 

where 

S d = ^l drdr'—±— (15) 
\l z Jn \r — r '| 

is a numerical constant which depends on the geometry. Thus repeating the summation of 
the diagrammatic series shown in Fig. |l](c) one obtains 

^=°> = T7(^' (16) 

where = S^e 2 ^, which is the usual random phase approximation (RPA) for the e-e 
interactions. Using the above result in Eq. (pi) results in 



/dN\ 



(17) 



\dfi/ \AJ \l + S d aU/nAJ \AJ \1 + {nL) d - 1 ) ' 

Therefore, the compressibility tends to its non-interacting value for kL 1 (i.e., large 
screening length, possible for example in very small semiconducting grains) and to (I/kL)^ 1 
of its non-interacting value for kL ^> 1 (metallic grains). 
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III. FLUCTUATIONS IN COMPRESSIBILITY 



In this section we shall calculate the fluctuations in the compressibility in the presence 
of e-e interactions. For the non-interacting electrons the fluctuations are expected follow 
the Wigner-Dyson statistics. As was noted in Ref. ||, in order to obtain the Wigner-Dyson 
statistics in the perturbative diagrammatic calculation one must insert a cut-off in energy 
(or temperature) of order of the single electron level spacing A. Explicitly the fluctuations 
can be calculated using the diagrams appearing in Fig. ^(a), which correspond to 

m q 

where u m = 2irmT is the Matsubara frequency and T is the temperature. D Wm (g) is the 
Fourier transform of the diffusion propagator which is the solution to the equation 



^) = ^EE^ m ®. (is) 



[OJr, 



DV 2 )2^ m (f,r = D5(f-f'), (19) 



with reflective boundary conditions Vf l T> LOrn (f, f ') = 0, where n is the normal to the sample 
edge. For a rectangular grain of dimensions L 3 , the solution is 

t> 1 V Ui=x,y,z cos (fc^) cos (Ay-) 

^u m \r,r ) - 2^ n„2 , ' y zyi > 

"n 1= -oo U 1 + UJ m 

where q 2 = k\^ + k 2 + k 2 g , (for a two dimensional system the z component drops out) and 
ki = irrii/L. After inserting the value of diffusion propagator given in Eq. ( pOl) one is left 
with the following summation: 

(* 2 7r) = T E E * + Dq 2 Y\ (21) 

\ C/i / m q 

The most significant contribution comes from the zero mode (q = 0) resulting in 

The fluctuations in the interacting case are represented by the diagrams shown in Fig. 
0(b) corresponding to 
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S 2 ^-)~ L ■ ^7777T7 ■ (24) 



5 2 |) = ^ 4 (o)EE^J?). (23) 

where x(0) is given by Eqs. (0,0)- Note that X appears in the power of four, which is the 
the result of the fact that when two RPA lines intersect, four lines appear. For the short 
range interactions we obtain 

1 

3/7/ " \l + a d U/QA J A*' 

Thus, the fluctuations in the compressibility are suppressed by the interactions as compared 
to the non-interacting value. Even the relative fluctuations (defined as the fluctuations in the 
compressibility divided by the averaged compressibility) are suppressed by the interactions. 
According to Eq. fllTI): 

: 2 dN\ //dN\ 2 ( 1 \ 2 

5 ^7//w/ ~{T^mK) ■ (25) 

A similar situation exists also for the Coulomb interactions. By using the value of x 
given in Eq. (|T6|) one obtains 

^9N\ ( 1 \ 4 1 



and (using Eq. (|T 



OAS > 2 ~ (ttoM • 



Thus, for kL <C 1 the relative fluctuations are the same as in the non-interacting case, while 
kL 3> 1 the relative fluctuations are suppressed by a factor of (KL)~ 2 ^ d ~ l > . It is important 
to note that for most metallic and semiconducting samples the above relation holds since k 
is of the order of the Fermi momentum, i.e., kL ~ 1 for a grain with ~ 1 electron. 



IV. CAPACITANCE 



From the experimental point of view, capacitance measurements is one of the most 
accessible methods to investigate the spectral properties of quantum dots and grains [||. 
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Therefore, it is interesting to connect the fluctuations in the compressibility to fluctuations 
in the capacitance for interacting disordered systems. The capacitance of a grain C may be 
related to its chemical potential in the following way: 



which is equivalent to 



when kL ^> 1 and thus one may write the average capacitance as 

<c> = (30) 

After inserting the calculated compressibility for the long-range Coulomb interaction (Eq. 
(|17|)) one obtains 

(C) = ±-L, (31) 

od 

which is the expected purely geometrical value of the capacitance. 

The fluctuations in the capacitance can be expressed via the fluctuations in the com- 
pressibility in the following way 
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{S' 2 C) = ( -n^- 1 A ) (5 2 ^), (32) 



,dN\ 



which after inserting Eq.(j27|) results in 



1 N 2(d-1) 

(5 2 C) = (C) 2 — . (33) 



Thus when kL ^> 1 the fluctuations in the capacitance are much smaller than the average 
capacitance. 



V. INVERSE COMPRESSIBILITY 



As mentioned briefly in the introduction, the inverse compressibility <9/i/ dN of a quantum 
dot can be deduced from measurements of the spacings between consecutive gate voltages 
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for which the conductance through the dot peaks. The spacings in the gate voltage are 
proportional to the difference between the chemical potential of N + 1 and iV electrons || 
denoted by A^, which may be related to the ground state energies in the following way: 

A 2 = fiN+i — = E N+ i — 2E N + Eff-i, (34) 

where En is the ground state energy of the dot populated by N electrons. In the continuum 
limit A2 = dfi/dN, i.e., the experiment actually measures the discrete limit of the inverse 
compressibility. 

The first step in connecting the results obtained from the compressibility calculations to 
the behavior of A 2 is to relate (dfi/dN) to (dN/dfi). With no loss of generality one can 
write 

<§M«f><)> 

and 

-2 



which if one assumes a well behaved distribution of the compressibility will result in 



(36) 



In a similar manner 



These relations are valid as long as (5 2 dN/dfi) < (dN/dfi) 2 . From Eqs. (pTT] , p~7D and (p4] , p6|) 



it can be seen that this condition is usually fulfilled since kL ^> 1. The above condition also 
holds in the non-interacting metallic regime for which the usual diagrammatic expansion is 



valid. [1T7L|TB|| Thus 
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for short range interactions and 

(|a Hl + (K£n w 1 + ^) A (40) 



\dN/ V v 1 J \ LA J 
for Coulomb interactions. The fluctuations in both cases are equal to 

i dp \ 



^)~A". (41) 
From Eqs. (|39l , f40l) it is clear that {dp/ ON) grows proportionally to the interaction strength 



U, while the fluctuations are independent of the interaction strength. 

This behavior implies that while the e-e interactions tend to shift the distribution of the 
inverse compressibility the width of the distribution will not change significantly fT9|| . 

The results for {dp/ ON) could also be anticipated from a simple assumption on the 
distribution of the electron density of the ground-state. The electrostatic energy needed to 
add an electron to a system of N electrons is given by 

eg t = f drdr' U(r,r ')p N (f) Pl (r '), (42) 

where p^ is the density of the N electrons already in the system, and p\ is the density of the 
additional electron. Assuming that both densities are uncorrelated, and that on the average 
{Pn) = N/fl and (pi) = 1/fi, the average electrostatic energy needed to add an electron is 

r N 
(eZ t )= j n dfdr>U(r,r>)-, (43) 

which for short-range interactions results in {sf n /} = Na d U/Q and for the Coulomb inter- 
actions {ef ni ) = NSde 2 L d ~~ 1 /Q = N{hcL) d ~ 1 A. Since under these assumptions {dp/dN) = 
( £ fntl ~ ^int 1 } + A> one immediately obtains Eqs. (|39|j40| ), which is the classical limit of the 
Coulomb blockade. |20| 



Also the fluctuations in the inverse compressibility may be deduced from similar assump- 
tions, resulting in 



<(*&) 8 > ~ <^> 2 , (44) 
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for both short and long range interactions. Thus, the fluctuations in the charging energy 
are (S 2 e^ lt ) = {(£f nt ) 2 ) — {^ n t) 2 ~ 0> This leads to the conclusion that there is no additional 
contribution to the inverse compressibility fluctuations beyond the fluctuations in A, in 



agreement with Eq. (f4T|) 



Therefore, we expect that the results presented in the previous sections will hold as long 
as the RPA assumptions hold, i.e., Vf 3> e 2 /h, where Vf is the Fermi velocity. 

VI. NUMERICAL CALCULATION OF THE INVERSE COMPRESSIBILITY 

In this section we shall numerically test the properties of the inverse compressibility. Of 
course, in a numerical calculation we can only compute the discrete form of dfi/dN, i.e, A2 
as a function of the strength of e-e interactions in the system. 

As a model system we chose a system of interacting electrons on a 2D cylinder of cir- 
cumference L x and height L y , which has been previously used in the study of the influence 
of e-e interactions on persistent currents. Ell] The model Hamiltonian is given by: 



H = e k,jal d a kd - Vj2( a l,j+i a k,j + a\ +lj a Kj + h.c) + H int (45) 

k,j k,j 

where a\ ^ is the fermionic creation operator, e k j is the energy of a site (k,j), which is chosen 
randomly between — W/2 and W/2 with uniform probability and V is a constant hopping 
matrix element. Hi nt is the interaction part of the Hamiltonian given by: 

H int = jJ2 a k,j a k,j0i,p a l,P ( 46 ) 
{k,j;l,p} 

for the short range interactions (where {. . .} denotes nearest-neighbor pair of sites) and 

for the Coulomb interaction, where b is the lattice constant. 

For a sample of M sites and N electrons, the number of eigenvectors spanning the many 
body Hilbert space is m = y^f). The many-body Hamiltonian may be represented by an 
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m x m matrix which is numerically diagonalized for different values of the e-e interaction. 
Here we consider a 4 x 3, 4 x 4 and 4x5 lattices with M = 12, 16, 20 sites correspondingly. 
For each value of M the average and fluctuations of A^ 2 are calculated. To obtain A^ 2 in 
each case Eg S for N = M/2 — 1,N = M/2 and N = M/2 + 1 are calculated. For the largest 
system considered (M = 20, N = 10) the calculation of the ground-state energy corresponds 
to diagonalizing a 184756 x 184756 matrix. We chose W = 8V for which this system is 
in the metallic regime pi) and average the results over 500 realizations for each value of 



interaction strength for the M = 12 and M = 16 cases and 200 realizations for M = 20. 

The results for (A^ 2 ) and ((<5 2 A^ 2 )) 1//2 are given in Fig. |3]a for the short-range inter- 
actions and in Fig. |3|b for the Coulomb interaction. It can be seen that in both cases (A^ 2 ) 
increases linearly for low values of U while (5 2 A^ 2 ) remains constant. This is in qualitative 



agreement with Eqs. (|39|- 41). In order to obtain also quantitative agreement one must take 



into account some finer details. For (A^ 2 ) one must consider the fact that the calculations 
were performed on a lattice. For short range interactions, one should replace a d /Q in Eq. 
fl3"5|) by Z(M)/(M — 1) where the average number of nearest neighbors, Z(M), depends 
on M due to the different ratio of sites close to the boundaries, where for M = 12, 16, 20, 
Z(M) = 3.333, 3.5, 3.6 correspondingly. Using the values of the single electron spacings 
A = 0.42V, 0.59V, 0.77V for M = 12, 16, 20 in Eq. (||) results in the curves plotted in Fig. 
1(a) for (Af /2 ). One can see a good fit up to U = V. For the long range interactions 
in our lattice system one should replace the integration in Eq. (|i"5D by a summation, i.e. 
S2 = (L/M 2 ) J2k,jjti, P l^fcj ~ ^,pl _1 > which results in S2 = 2.35,2.44,2.48 corresponding to 
M = 12, 16, 20. After incorporating S2 and A in Eq. (flTf ) one obtains the curves plotted in 
Fig. |b. An excellent fit is seen up to U = 3V. The exact value of (5 2 Af /2 ) = (4/vr - 1)A 2 
is deduced from RMT and plotted in Fig. [| It can be seen that this prediction is also 
confirmed for weak interactions p2| . 

It is possible to observe what is the influence of the e-e interactions on the full distribution 
of A 2 . In Fig. ^ the distribution for the non-interacting case (U = 0) is compared with the 
distribution for U = 2V in the long range interaction case. As can be seen in Fig. [|b, for this 
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value of interaction strength no significant change in the second moment of the distribution 
(<5 2 A^ 2 ) is expected, and indeed the width of the distributions is almost identical. On 
the other hand, higher moments seem to be influenced by the e-e interactions resulting in 
changes in the tails of the distribution. The effect of interactions on the higher moments of 
the A^ 2 merits further studies. 

As has been mentioned in Ref. ||, and discussed in the previous section, we expect the 
analytical results based on the RPA approximation to hold as long as no correlations develop 
in the electron density. Following Ref. || we define a two point correlation function: 

(jO-^ _ ^k,j>l,pC(r k j — n.p^lrfc.j-n.pl.r 
^2k,j>l,p ^|r fej -n iP |,r 

where 

r(if * s _ / [gl^fcj ~ (^kM^p^P ~ ("Lai*)] \ (AQ , 

\ ^k,j a k,j){al p ai,p) I 

In Fig. |] we present the correlation C(r = y2b) (which corresponds to a pair of diagonal 
sites) for the long-range interactions as function of the interaction strength. Under the RPA 
assumptions we expect C(r — V2b) — ► 0. As can be seen there is some correlation as result 
of the boundaries even without e-e interactions. As the interaction strength increases there 
is no dramatic change in C(r = y/2b) up to U = 2.5V and then one sees a strong deviation. 
A similar behavior can be seen for the short range interactions where the deviation appear 
at U = V. This agrees well with the values of interaction for which the numerical results 
depart from the RPA analytical calculations. 

One can roughly estimate the strength of interaction U for which these correlations should 
appear. As mentioned in the previous section, we expect the RPA approximation to hold 
while e 2 /vf < 1, which for long range interactions corresponds to r s = yj V / 2(U / AV) < 1 



and for short range interactions to r s = ^ir/2(ZU/AV) < 1. This is surprisingly close to 
the values for which deviations from RPA theory are seen numerically. 
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VII. DISCUSSION 



The main conclusion that follows of the previous sections is that the RPA approximation 
fails to explain the large fluctuations in the inverse compressibility seen in the experiment 
0-0. Under the assumptions of an RPA treatment of the e-e interactions the fluctuations 
are proportional to the single electron mean level spacing, therefore independent of the 
strength of the e-e interactions. In contrast, experiments show that the fluctuations are 
substantially larger than the mean level spacing, and seem to be about 15 percent of the 
average inverse compressibility. 

Actually, one could have anticipated the failure of the RPA calculation for the experi- 
mental realizations, since their densities are too low for the RPA approximation to remain 
valid. For a 2DEG quantum dot one may rewrite the condition for which the RPA fails 
r s = e 2 /vf < 1 as n < l/ira 2 B , where n is the electron density in the dot and a# is 
the Bohr radius, which in GaAs is ~ lOOA. Thus one expects RPA to fail at densities 
?i<3x 10~ 11 cm -2 . In all of the recent experiments for which large fluctuations were ob- 
served |3|-|5| the 2DEG density is about 3 x 10~ n cm -2 , and the density in the dot is probably 
even lower. 

The large fluctuations in the inverse compressibility are caused by spatial fluctuations 
in the electron density which are not taken into account in the RPA approximation [J5J. It 
is important to note that although the charge distribution becomes inhomogeneous when 
the electron density is still much higher than the one critical for Wigner crystallization. 
Thus the experiments are apparently in an intermediate range of densities (or intermediate 
e-e interaction strength). In this regime RPA approximation does not hold anymore and 
some density correlations do appear, however the systems are still very far from the Wigner 
crystallization. 

Note that in the case of large dimensionless conductance g ^> 1 that was considered 
in this paper we found no influence of the disorder on the structure of Coulomb blockade 
peaks. For different physical reasons fluctuations of the order of the inverse compressibility 
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may also appear in the localized regime. |23| 



Finally, it is also interesting to compare the fluctuations in the compressibility to the 
fluctuations in the polarizability. ]15[ Although the calculation of both quantities show many 
similar features, there is one crucial difference. The main contribution to the fluctuations 
in the compressibility comes from the zero-mode, which corresponds in the non-interacting 
case to the single electron level fluctuations. In the polarizability, there is no contribution 
from the zero mode, and the fluctuations stem from the other modes, which correspond to 
density fluctuations. This results in the fact that the relative polarization fluctuations are 
smaller than the relative compressibility fluctuations by a factor 1/g 2 . 

ACKNOWLEDGMENTS 

We would like to thank A. Auerbach, M. Field, Z. Ovadyahu, B. Shklovskii, U. Sivan 
and D. A. Wharam for useful discussions. R.B. would like to thank the US-Israel Binational 
Science Foundation for financial support. 



16 



REFERENCES 



[1] For recent reviews see: M. A. Kastner, Rev. Mod. Phys. 64, 849 (1992); U. Meirav and 
E. B. Foxman, Semicond. Sci. Technol. 10, 255 (1995). 

[2] R. C. Ashoori, Nature 379, 413 (1996). 

[3] A. M. Chang, H. U. Baranger, L. N. Pfeiffer, K. W. West and T. Y. Chang, Phys. Rev. 
Lett. 76, 1695 (1996). 

[4] J. A. Folk, S. R. Patel, S. F. Godijn, A. G. Huibers, S. M. Cronenwett, C. M. Marcus, 
K. Campman and G. Gossard, Phys. Rev. Lett. 76, 1699 (1996); C. M. Marcus, J. 
A. Folk, S. R. Patel, S. M. Cronenwett, A. G. Huibers, K. Campman and G. Gossard 
(preprint). 

[5] U. Sivan, R. Berkovits, Y. Aloni, O. Prus, A. Auerbach and G. Ben-Yoseph, Phys. Rev. 
Lett. 77, 1123 (1996). 

[6] V. Chandrasekhar, Z. Ovadyahu and R. A. Webb, Phys. Rev. Lett. 67, 2862 (1991). 

[7] L. P. Gorkov and G. M. Eliashberg Zh. Eksp. Teor. Fiz. 48, 1407 (1965) [Sov. Phys. 
JETP 21, 940 (1965)]. 

[8] B. L. Altshuler, and B. I. Shklovskii, Zh. Eksp. & Teor. Fiz. 91, 220 (1986) [Sov. Phys. 
JETP 64,127 (1986)]. 

[9] M. L. Mehta, Random matrices (Academic Press, San-Diego, 1991). 

[10] K. B. Efetov, Adv. Phys.32, 53 (1983). 

[11] J. J. M. Verbaarschot, H. A. Weidenmuller, and M. R. Zirnbauer, Phys. Repts. 129,367 
(1985). 

[12] R. Berkovits, Europhys. Lett. 25, 681 (1994); R. Berkovits and Y. Avishai, J. Phys. 
Condens. Matter 8, 391 (1996). 



17 



[13] S. Strassler, M.J. Rice and P. Wydner, Phys. Rev. B 6, 2575 (1972). 

[14] M. J. Rice, W. R. Schneider, and S. Strassler, Phys. Rev. B 8, 474 (1973). 

[15] R. Berkovits and B. L. Altshuler, Europhys. Lett. 19, 115 (1992); Phys. Rev. B 46, 
12526 (1992). 

[16] B. L. Altshuler and A. G. Aronov, in Electron- electron interactions in disordered systems 
Ed. A. L. Efros and M. Pollak (North-Holland, 1985) 

[17] B. L. Altshuler, Y.Gefen and Y. Imry Phys. Rev. Lett. 66, 88 (1991); A. Schmidt (1991). 

[18] A. Kamenev, B. Reulet, H. Bouchiat and Y.Gefen, Europhys. Lett. 28, 391 (1994). 

[19] On the other hand, since we did not calculate higher moments, it is possible that the 
tail of the distribution will change as function of the interaction strength. 

[20] D. V. Averin and K. K. Likhraev in Mesoscopic phenomena in solids Ed. B. L. Altshuler, 
P. A. Lee and R. A. Webb (North-Holland, 1991) 

[21] R. Berkovits and Y. Avishai, Europhys. Lett. 29, 475 (1995); Phys. Rev. Lett. 76, 291 
(1996). 

[22] Some deviations occur for M = 20, probably because for this larger lattice size the 
system is closer to the Anderson transition. 

[23] A. L. Efros and B. Shklovskii, in Electron- electron interactions in disordered systems 
Ed. A. L. Efros and M. Pollak (North-Holland, 1985); F. Pikus and B. Shklovskii (APS 
March meeting abstract, 1996) 



18 



FIGURES 



(a) 



(b) 




(c) 

V ) 



FIG. 1. Feynman diagrams representing the average (a) non- interacting and (b) interacting 
compressibility., where the wavy lines correspond to the electron-electron interaction given in dia- 
gram (c). 
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FIG. 2. (a) Feynman diagrams representing the fluctuations in the compressibility of a sample. 
Diagram (a) corresponds to the non-interacting case, while diagram (b) to the interacting case. 
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FIG. 3. (Af /2 ) and (<5 2 Af /2 ) for different lattice sizes (3x4,4x4 and 5 x 4) as function of 
interaction strength for (a) short range interactions (b) Coulomb interactions. The lines represent 
the theoretical predictions given in the text. The full line corresponds to a 3 x 4 lattice, the dotted 
line to 4 x 4 and the dashed line to 5 x 4. 
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FIG. 4. The full distribution P(Aj /A), for a 4 x 4 lattice. The full line corresponds the GOE 
distribution P(x) = (ttx/2) exp(— ttx 2 /4), where x = A^ 2 /A, and the dashed line corresponds to 
the GOE distribution shifted by the average charging energy for U = 2V, i.e., by A^ 2 +A— (A^ 2 ) . 
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FIG. 5. The two point correlation C(r = 2 l l 2 b), for long range interactions, as a function of 
the interaction strength. 
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